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ABSTRACT: The double-homodyne and the heterodyne detection schemes for phase shifts between 
two synchronous modes of the electromagnetic field are analyzed in the framework of quantum 
estimation theory. The probability operator- valued measures (POM 5 s) of the detectors are evaluated 
and compared with the ideal one in the limit of strong local reference oscillator. The present 
operational approach leads to a reasonable definition of phase measurement, whose sensitivity is 
actually related to the output r.m.s. noise of the photodetector. We emphasize that the simple- 
homodyne scheme does not correspond to a proper phase-shift measurements as it is just a zero-point 
detector. The sensitivity of all detection schemes are optimized at fixed energy with respect to the 
input state of radiation. It is shown that the optimal sensitivity can be actually achieved using suited 
squeezed states. 


1 Introduction 

Weak forces on macroscopic bodies in interferometric arrangements or, more generally, minute variations 
of environmental parameters in optical fibers are detected through the induced changes in the optical 
paths of the light beams. The detection of the induced phase-shift represents one of the most sensitive 
measurement on radiation in order to monitor such small perturbations. The back-action effect on the 
measured parameter due to the radiation pressure imposes limitations on the radiation intensity, and 
improvements of the sensitivity can only be achieved by suited preparation of the input quantum state. 

In this paper a narrowband analysis of some relevant phase detection schemes is presented (a multimode 
wideband analysis can be found in [1]). Classically a phase- shift measurement in an interference experiment 
can be directly related to the polar angle between two quadratures of one field-mode, which in turn are 
given by two output photocurrents. Quantum mechanically the quadratures of the field are noncommuting 
observables and their relative polar angle cannot be interpreted as a selfadjoint operator, as in the early 
Dirac’s heuristic approach [2]. The oscillator phase is not an observable in usual sense, and the problem 
of identifying its quantum dynamical counterpart has provoked many discussions in the literature (see, for 
example, [3, 4, 5] and references therein). Among the numerous attempts the limiting procedure of Pegg 
and Barnett [3] has become the most popular technique, because it allows the evaluation of expected values 
with very simple and reliable rules. However, despite its simplicity and effectiveness as a mathematical 
tool, this approach has no obvious physical interpretation, and leaves most of the conceptual problems 
on phase detection still open. The actual problem of a phase measurement description does not concerns 
with a definition of a selfadjoint operator, but with an operative recipe to evaluate the phase statistics in 
a real measurement, starting from the knowledge of the density matrix of the input radiation. On these 
lines the most appropriate approach to the phase detection of the field is the quantum estimation theory 
of [4]. Even though it easy to show that this method is equivalent in the end to the Pegg and Barnett 
procedure (see [8] for more details), nevertheless it provides a physically meaningful scheme for the phase 
measurement where all conceptual problems disappear. Despite it has long been recognised as the most 
natural framework for analyzing any kind of quantum detection, the quantum estimation theory has not 
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gained the necessary popularity yet, perhaps due to the fact that its main ingredient— the probability 
operator measure (POM) — is generally a nonorthogonal spectral decomposition, and thus appears to be 
in conflict with the conventional dictum of quantum mechanics that only “observables ” i.e. orthogonal 
POM’s— can be measured. This point has been well clarified in some papers (see for example [8]), where 
it is shown that nonorthogonal POM’s correspond to actual observables on a larger Hilbert space which 
includes also modes pertaining the measuring apparatus (all together referred to as “probe”). It is clear 
that this assertion provides the proper quantum setting for the operational point of view of [11], where the 
dependence of the measured operator on the detection scheme just corresponds to the involvement of the 
probe variables in the measuring process, involvement which becomes unavoidable when the phase of the 
field is detected. 

Any quantum measurement needs a classical final stage and for measurement on radiation this is essentially 
a photodetection. Moreover a proper measurement of the quantum phase has to be related to the detection 
of a quantity which itself is a phase, i.e. is defined on the unit circle. In this sense we distinguish between 
two main different classes: the genuine phase detection schemes and the measurements of a single phase- 
dependent observable. In the former class, the phase-shift of the field is related to the polar angle between 
two measured photocurrents which, in turn, correspond to a couple of two conjugated quadratures of 
the field. Such scheme is the only viable one for phase detection, and corresponds equivalently to either 
heterodyning or double-homodyning the field. This also clarifies the subtle nature of the phase itself 
which, despite being a single real parameter, nonetheless requires a joint measures of two conjugated 
operators. In contrast, in the second class of measurements, only a single observable is detected— typically, 
when homodyning a single-quadrature of the field. Here we want emphasize that a single- quadrature 
measurement cannot be used to infer the value of the phase, because the knowledge of a quadrature would 
require an additional measurement on the field— essentially its intensity— which unavoidably would destroy 
the information on phase. Thus, the single-homodyning scheme can be used only as a zero-phase monitoring 
technique, which, however, is the essential of a typical interferometric measurement. In order to stress the 
operational nature of POM approach here we also present, as an example, a numerical simulation of a real 
experiment, reproducing the classical output photocurrents due to a particular quantum state, and then 
evaluating the phase statistics as the polar angle distribution. 

After selecting a measurement procedure — either ideal or feasible — the statistics of the detected phase can 
be further improved at fixed total energy by looking for optimal states of the field. We show that the r.m.s. 
phase sensitivity versus the average photon number n is bounded by the ideal limit A<£ ~ n , whereas 
for the feasible schemes the bound is A <j> ~ n _2/3 , in between the shot noise level A <j> ~ n 1/2 and the 
ideal bound. The latter can actually be achieved by single-homodyning suitable squeezed states, but only 
in the neighborhood of a fixed zero-phase working point. The requirement of detecting the whole phase 
probability distribution makes the proper phase measurement less sensitive than in the case of a zero-point 
detection and the state achieving the two bounds are dramatically different: they are weakly squeezed 
(about 2% of squeezing photons) for the double- quadrature measurement, whereas they become strongly 
squeezed (50%) for one-quadrature detection. 

Sect.2 is devoted to the theory of phase measurements, with a detailed analysis of the various schemes. 
Subsect. 2.1 is a brief review of the quantum estimation approach. Subsect. 2.2 presents some remarks 
and criticisms about the different definitions of sensitivity. Detection of the phase through simultaneous 
measurement of two quadratures of the field are discussed in Subsects. 2.3 and 2.4, which are dedicated 
to double-homodyne and heterodyne detection (it is shown that the two schemes are fully equivalent). 
Subsect. 2.5 examine s the measurement of the quantum sine and cosine of the phase, with a comparison 
between the present quantum estimation approach and earlier treatments [7] .In Subsect. 2.6 we analyze the 
homodyne detection scheme. In Sect.3 the optimal states of the field are given, which maximize sensitivity 
for all the schemes of Sect.2, also indicating how to actually achieve such states. 
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2 Quantum measurement of the phase 


At presently it is fairly universally accepted that the quantum description of a phase measurement on 
a single mode of the e. m. field cannot be approached by means of the usual concept of observable. 
As a matter of fact, even though selfadjoint operators can be defined on the Fock space, none of them 
can appropriately describe the actual statistics of phase measurements. The quantum estimation theory 
of Helstrom [4] provides the most general — nonetheless operational — framework to analyze any kind of 
measurement in the quantum context, and, in particular, the detection of a phase differences between 
synchronous oscillators. The main ingredient of such theory is represented by the mathematical concept 
of probability-operator- valued measure (POM) on the Hilbert space Tis of the system, which extends the 
conventional description by selfadjoint operators. Using a notation which is familiar to physicists — even 
though not strictly legitimate from the mathematical point of view [8] — given a set of (generally complex) 
parameters z to be measured, a POM dX(z) is a self-adjoint measure with the following properties 

dX(z) > 0 , f dX(z) = i , (2.1) 

Jz 

where Z denotes the space of the parameters z. Eqs.(2.1) assure that using a POM one has an operational 
recipe which univocally relates the density-matrix state p of the system to a probability distributions of 
the parameters z corresponding to a particular experimental setup. In formulas one has 

dP{ z) = tr{pdX(z)} . (2.2) 

Also a set of selfadjoint operators A can be defined 

A = [ zdX(z) , (2.3) 

j z 

and, more generally, operator functions /(A) as 

/( A ) = f /(z)^(z) • (2.4) 

When the POM dX{z) = |z)(z| dz is given in terms of orthogonal |z)’s, then it corresponds to the customary 
measure of the commuting observables A, whose corresponding selfadjoint operators obey the function 
calculus 

fW = /(A) . (2.5) 

On the contrary, the relation (2.5) no longer holds true for a generic nonorthogonal POM. As a consequence, 
the selfadjoint operators A only provide the expected values of the parameters z, whereas the higher 
moments of the probability distributions differ from the corresponding moments of the operators A 
themselves, and can be evaluated only through Eq.(2.2). One should notice that, despite the POM’s 
generally describe measurements that do not correspond to observables in the usual sense, nonetheless 
there is no conflict with the basic assertion of quantum mechanics that only observables can be measured. 
In fact, the Naimark theorem assures that every POM can be obtained as a partial trace of a customary 
projection- valued measure on a larger Hilbert space [8] which itself represents the original system interacting 
with an appropriate measuring apparatus. Upon denoting by |^( z )) G 7ts 0 Tip a complete orthonormal 
set of eigenvectors of commuting selfadjoint operators acting on the enlarged space — including also the 
apparatus (probe) space Tip — the probability distribution 

dP{z) = tr(ps® pp\ip{z))(rp(z)\) = tr s [pstr P (p P \^(z))(tp(z)\)] , (2.6) 
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corresponds to the POM on the Hilbert space 7is of the system only 

d\(z) = trp (p/»M*)XlH*)l) = (V r (*)I^HV , (*)) • ( 2 - 7 ) 

Notice that the above extensions of the system Hilbert space— hence, physically, the experimental 
realisation of the POM — is not necessarily unique. 

In Subsect. 2.1 we review the quantum estimation theory of the phase which leads to the optimal POM 
representing the most accurate measurement. This, however, is only an ideal limit as no viable schemes 
impl ementing such POM has been envisaged yet. Therefore, in Subsects. 2. 3 and 2.4 experimentally 
achievable detection schemes which correspond to a sub-optimal POM (double- homodyne and heterodyne) 
are analysed in detail, whereas in Subsect .2.6 the customary homodyne detection is considered. The latter 
exits from the present quantum estimation approach, however it is in order, due to the relevance of this 
scheme in any interferometric setup. 

2.1 Canonical Measurement 

The quantum estimation theory analyzes the possible strategies for estimating a parameter on the basis 
of an error-cost-function: the optimum POM is the one which minimizes the total average cost. For a 
maximum-likelih ood criterion, the optimum POM for the phase is 

dm = . ( 2 - 8 ) 

\e'+) being the Susskind-Glogower phase states [7] 

|e*) = £ e ‘"* I”) • ( 2 - 9 ) 

n=0 

It is worth noticing that in the present case the maximum likelihood criterion is equivalent to the Gaussian- 
cost-function one for the bounded case (sin 2 - cost) [4]. Such generality of the optimal POM dji(<f>) justifies 
the term Canonical Measurement here adopted for this POM approach. Also notice that the Pegg and 
Barnett [3] approach is totally equivalent to the present one as regards evaluations of statistics. However, 
there is no physical interpretation for the mathematical tricks on which their method relies. Some examples 
of commuting pairs of self-adjoint operators achieving the optimum POM (2.8) on a system-probe Hilbert 
space have been proposed in [9] and in [5]: however, no viable method for experimentally implementing a 
corresponding setup has been devised yet, and hence the POM (2.8) only represents an ideal limit. 
Corresponding to the optimal POM (2.8) one defines the selfadjoint phase operator 

f = f * dm = -i£ (-r™^ MW ( 21 °) 

" w n^m 

and the squared operator 

? = /’ dm = J + 2 £ (-r'An MMI • (2-11) 

J * n^m * 

Notice that, as announced, 

<t> 2 ^ <t? , (2.12) 

and more generally f(<f>) ^ /(<£). Such fail of the operator function calculus also holds true for the 
selfadjoint operators defined through the experimentally feasible non-optimal POM’s. The fact that there 
is no orthogonal optimum POM for the phase, physically corresponds to the impossibility of defining the 
measurement of the phase independently on the apparatus. This assertion clarifies and formalizes the 
operational nature of the phase detection which has been pointed out by Mandel et al.in [11]. 


132 



2.2 Phase Sensitivity 

Usually the sensitivity of a measurement of a parameter say x e R is assumed equivalent to the r.m.s. of 
the experimental probability distribution dP(x ) , namely 

Ax 2 = f R dP{x)x 2 - (f R dP(x)x ) . (2.13) 

On the other hand, the phase variable <f> is defined in the bounded domain [-*■,*•) with 2x-periodicity: 
this peculiar property of the phase has lead many authors to the conclusion that the r.m.s. of the phase is 
not the appropriate quantity to be considered as an evaluation of the phase sensitivity, claiming that it is 
not invariant under phase variable translation (j> — » <f> + x . Thus, different definitions for the sensitivity 
have been adopted, which would be equivalent for an unbounded Gaussian- distributed variable. Here, 
after critical reviewing such quantities, we show that an operational definition of a phase measurement 
procedure leads unequivocally to adopt the r.m.s. itself as the correct sensitivity parameter. 

a. Phase dispersion D[ 12, 13]. 

Dispersion D is defined as follows 


D = (1 — (cos <f>) 2 - (sin 0} 2 )2 


1 - 


]^ C »i+l C n > 
n— 0 


(2.14) 


where c n are the coefficient of the number representation of the state and the sine and cosine operators are 
defined according to Eq.(2.4) as follows 


cos 


<t> = [ cos <f> dfi{4>) , 
J — IT 

~<t> = J sin <f> dfi((f>) , 


(2.15) 

(2.16) 


and coincide with the sine and cosine operators of Susskind and Glogower [7]. The definition (2.14) 
follows from elementary error-propagation calculus, the phase <f> being regarded as a function of the two 
"independent variables” sin <j> and cos <j> as follows 


<f> = — tln(cos^ -f isin<£) . (2.17) 

In Eq.(2.17) the correct logarithm branch is selected in order to obtain the desired domain for <f>. A part 
the m inor point that Eq.(2.14) would lead to dispersion D — 1 for constant distributions — instead of 
~ * ma in criticism is that sin <f> and cos <f> cannot be considered as independent variables, 

because they correspond to a noncommuting pair of operators which are jointly measured when detecting 


b. Reciprocal peak likelihood £^[5]. 

The peak likelihood P(<f> \<f>) is the maximum height of the probability distribution. Its inverse, namely 


6<t> = 


1 

P(^) 



(2.18) 


has been introduced in [5] as a measure of the width of the distribution, coherently with the mtuxxmum - 
likelihood strategy used in the quantum estimation theory. Here, the following criticisms are in order: i) S<f> 
is a local criterion, namely it checks only one point of the distribution, whereas there is no control on the 
global behaviour as, for example, on the eventual occurrence of high tails. The most degenerate situation 
occurs when the tails are so high that the distribution itself converges to the P{<f>) = 1/2t apart from 
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one point with infinite probability density and zero integral [14], thus leading to vanishing &<j> instead of 
66 = x 2 /3; ii) the coherence of this sensitivity definition with the maximum-likelihood strategy [5] cannot 
be considered as a valid argument, in view of the aforementioned equivalence between the likelihood strategy 
and the (quasi)- Gaussian one; iii) recent numerical results [6] have shown that the simulated sensitivity 

does actually not correspond to 6<t>- 

Given f appittos (or an ideal detector) one has a corresponding POM and, in turn, a probability 

distribution dP{6) according to Eq.(2.2). Such probability has a r.m.s. error (2.13) given by 


(A**) = ( 6 2 ) ~ ( 4 >? 


(2.19) 


Here (. . .) denotes the ensemble quantum average on the system space H s , and . the °P er * tor * * “ d ** 
depend on the considered POM (for the optimum POM they ^are given m Eqs.(2.10,2.11)). Notice that 

there is no ambiguity in choosing between the two operators 6 2 * because Eq.(2.19) 

from the probability (2.2). For a random-phase state— namely a constant probability distribution— one 

tStoofL»bmc« under ph»e drift., here we .tre.. that thi, eetudly i. not . problem, 
to fact, the only concern i» the corre.pondence between experimental and theoretical quantitie^and the 
circular topology of the phase arises at both experimental and theoretical levels in the same way. Whae™* 
procedure is considered for measuring the phase, the information on it has always to be ^ erred ^ 
joint sine-cosine measurement, and hence the experimental equipment itself has to be tuned on a selected 
2x-window. Once the domain is fixed, the experimental noise is, by definition, the r.m.s. noise on su 
domain. Therefore, different choices of the 2x-window actually lead to different experimented amount, of 
noise, and also theoretically the r.m.s. noise has to be evaluated on the chosen domain (hereafter we will 

always use the [— x, x) window). 

2.3 Double-homodyne Detection 

The double-balanced-homodyne [11] (DBH) detection provides a way for simultaneously measuring a couple 
of field-quadratures for one mode of e. m. field. The schematic diagram of the experimental set-up 
reported in Fig. 1. There are four 50-50 beam splitters and four identical photocounters, and a x /2 Phase 
shifter is inserted in one arm. The mode supporting the phase is a, whereas a stable reference for 
phase is provided by a local oscillator (LO) which is synchronous with a and is prepared m a highly exci 

The^DBH^eme can also perform a phase measurement, however with a probability , distribution ^ which 
does not correspond to the ideal case due to unavoidable addition of "instrumental noise. The DBH 
phase distribution is obtained through the following procedure. Each experimental event consists o 
simultaneous detection of the two difference photocurrents Ji = *e - »s and J 2 = n 4 - n 3 whidi trace 
two field-quadratures. Each event thus corresponds to a point plotted in the complex plane of the field 
amplitude. The phase relative to the event is nothing but the polar angle of the point itself An 
histogram of the phase distributions is thus obtained upon dividing the plane into nidi _( infinitesimal ) 
angular bins of equal width 66, bom -x to x, then counting the number of points which faU into each bin. 
In formulas, one has the statistical frequency P n for the n-th bin 0 n = [-x + n66 , -x + (» + 1)66) 


p n = -L of events with I\ = p cos 6 » h = P sin 6 , ^ € On} > 


( 2 . 20 ) 


where p = Jit + and N is the total number of experimental points. 

In Fig .2, as an example, a computer simulation of the above experimental procedure is iliustrated fora 
squeezed state with equal number (n) = 10 of signal and squeezing photons. The experimental histogr 
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Figure 1: Outline of scheme of a double-homodyne detectors 


(10 4 events) is compared with the theoretical results from the POM for the DBH detection. This can 
be obtained as follows. The difference photocurrents I\ and I 2 are commuting operators with factorized 
probability P(/i,/ 2 ) = P(Ji)P(/ 2 ). Introducing the reduced current T = // |z| for each homodyne detector, 
one has the probability distribution in terms of the Fourier-transform of the generating function for the 
moments (e*^) 


pm = /' w , ^ 

J — IT * ^ 


'*"■1*1 d\ 

n 1*1 

The phase distribution is the marginal probability integrated over the modulus p 


P{4>) = [ pdp P\ (p cos 4>)P 2 {p sin <f>) . 

Jo 

Using Eq.(2.21) one has 

P(4)= f°° pdp r p r tT ^ s ® P p 

JO J~t\z\ 2t /—it | x| 

ps being the density matrix of the mode a (the system) and 

pp = |0)(0| ® |0)(0| ® \z){z\ 


( 2 . 21 ) 


( 2 . 22 ) 


(2.23) 


(2.24) 


the density matrix of the probe. From Eqs.(2.6,2.7) one can see that the "experimental” POM is obtained 
upon tracing over the probe Hilbert space Jip y thus obtaining the operator which acts solely on the system 
space Hs 


dp D {(j>) = d<t> r pdp r H ^ r w ^ t TP {\ s ® p P 

Jo J-k\z\ 2x y_ K | z | 2 jt 1 


(2.25) 
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i t 

Figure 2: Computer simulation of the DBH experimental procedure for a squeezed state with equal number 
(n) = 10 of signal and squeezing photons. The experimental histogram (10 4 events) is compared with the 
theoretical results from the POM in Eq.(2.26) 


Using the coherent state resolution of the identity, the following closed formula is obtained (the detailed 
derivation is reported in the appendix) 


<M*> = * £ |„>H , 


(2.26) 


where T(x) is the Eulero’s gamma function. 

The POM in Eq.(2.26) for the DBH detector corresponds to an effective measured phase operator which 
is given by 

* f . .. / / \n-m 1 I \ / I (2.27) 


<t>D = J <t> d(l D {<t>) = -* Yj (“)" m — 


m -v/n!w»! 


I")M , 


and the squared one 


« - 1 * *** = t + 2 £ (-r vV*%^ ww • 


needed for evaluation of the instrumental sensitivity (A<^). For any state of the mode a one can simply 
verify that 

(A? d ) > (A? 2 ) , ( 2 - 29 ) 

namely the DBH scheme adds extrinsic instrumental noise, as it does not implement the optimal canonical 
measurement of the phase. However, we stress again that the DBH detection is the best available method 
for detec ting the phase. In Fig.3 a comparison between the canonical (ideal) and DBH (feasible) phase 
probability distributions is given for the same state of the computer simulation in Fig.2, showing that the 
former is sharper and higher than the latter. 
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Figure 3: Comparison between the ideal and double-homodyne phase probability distributions for the same 
squeezed state of Fig.2 

2.4 Heterodyne Detection 

The first proposed method to perform simultaneous measurements of two field- quadratures was the 
heterodyne detection. Here we synthetically analyze this scheme, only in order to make a connection 
with the double- homo dyne detector and show that that the two apparatus are completely equivalent from 
the point of view of the measured physical quantities. The input field E IN impinges into a beam splitter 
and has nonzero photon number only at the frequency u>o + u>j/r. The local oscillator works at the different 
frequency u>o, and the output the photocurrent Iout Is measured at the intermediate frequency w if . The 
measured photocurrent is given by 

IoUT{t) ~ ■®oi/t(0-®oi/t(0 » (2.30) 

where E ± denote the usual positive and negative frequency components of the field. The component at 
frequency u >jp is given by 

Iout{uif) - J dt iouT{t)e' u,Fl 

= J <L)E out (u> + ui > if)Eq UT (u) . (2.31) 

For a nearly transparent beam splitter, and in the limit of strong LO in the coherent state \z ) one can 
define the reduced complex current y 

y = lim 7 ~ 1 Iout{vif) , (2.32) 

r? -> 1, |z| — oo 

7 = cost. 
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where 7 = |zl\Aj(l - v)- la this limit the expression of y is given by 

y = |z| _1 (aj6j + a,6|) + vanishing terms , (2.33) 

where the subscript s, l and * refer to the signal, LO and image component of the field respectively, a are 
■ign^l modes, b the LO modes, and the vanishing terms denote operators which do not give contributions 
in the strong LO limit. In the double-homodyne in the same limit the role of the complex current (2.33) 

is played by 

i = ii + & 2 = |z| -1 (02a! + *o<4) » ( 2 * 34 ) 

where subscript 1 refers to the input signal and subscript 2 to the local oscillator, whereas to i» the vacuum 
mode at the unused port of the beam splitter which contains the input signal. The fully equivalence 
between heterodyne and double-homodyne is apparent when comparing Eq.(2.33) and Eq.(2.34). As in 
the double-homodyne case, now the real and imaginary parts of the current trace the two conjugated 
quadratures a* and a^ x/2 of the signal mode. In [16] the POM of the heterodyne detector has been 
derived in a different context, leading to the same result obtained for the double-homodyne in Subsect.2.3. 
Wc notice that the actual sources of extrinsic added noise are the vacuum modes a, for the heterodyne 
detector and 60 for the double-homodyne: the other vacuum modes are totally irrelevant in the limit of 

strong LO. 


2.5 Measurement of the phase quadratures 

The POM approach naturally leads to well defined operator functions of the phase which obey the 
trigonometric calculus at the operator level, and, hence, also at the level of expectation values. In particular, 
the sine and cosine operators are defined as in (2.16). Such definitions coincide, in the case of optimum 
POM, with the sine and cosine operators s and c introduced by Susskind and Glogower [7] 



c=l(e + e+) , 


(2.35) 


where i± denote the raising and lowering operators e+|n) = |n + 1), e_ = (e+)t. Notice, however, that 
ttiU equivalence between operators fails for higher powers, namely 

c" ± cos" ^ , i" f sin"^ , for n > 1 (2.36) 


Here some remarks are in order, regarding relevant differences between a conventional measurement of a 
single phase-quadrature— say the cosine c— and a joint measurement of both sine-cosine quadratures which 
have been analyzed in previous Subsections. A single phase-quadrature measurement leads to violation of 
the trigonometric calculus for expectation values. In fact, for a general density matrix state p one has that 


Tr[p(c 2 + i 2 )] = l-^(0|p|0), 


(2.37) 


whereas for a joint measurement one obtains 

Tr[p(sin 2 ^ + cos 2 ^)] = 1 . (2.38) 

We stress again that, however, the linear operators coincide in the two cases, and thus one gets the 
same average values. However, the probability distribution of the outcomes from single phase-quadrature 
measurement exhibits unphysical features for nonclassical states, whereas the probability distribution from 
the joint measurement does not. In the single- quadrature measurement one has 

P{c) = tr{p|c)(c|} , ( 2 - 39 ) 
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where the eigenstates of c are given by [7, 15] 


|c> = ^~(1 - c 2 ) ^ sin[(n + 1) arccos c]|n) 

n=0 

On the other hand, the Radon- Nikodym derivative of the joint measurement POM’s leads to 
P(c) = Tr 


P^yP-^] = ^(1 - c 2 ) > ^(m|^|n) exp [*(n - m) arccos c] , 

n,m 


d<f> dcj X 

for the opti mum POM case, whereas for double-homo dyning one obtains 


(2.40) 


(2.41) 


P(c) - Tr = 1(1 _ c 2 ) > X)(m|p|n) r ^ 2 + ^ exp [t(n - m) arccos c] . (2.42) 

^ * n m V Tl*77i. 


The differences between the single- quadrature and double- quadrature probabilities become striking for 
isotropic states, as, for example, the vacuum or a general number state. In this case the above distribution 
should be compared with the Radon-Nikodym derivative of the constant distribution 


P{c) = 


1 1 
tVI-c 2 ' 


(2.43) 


which is a concave function and has poles at the c = ±1 stationary points of the cosine. The probabilities 
(2.41) and (2.42) coincide with (2.43) for number states, whereas the probability (2.39) has the opposite 
curvature for the vacuum state, and oscillates fastly around the function (2.43) for nonvacuum n umb er 
states. These undesired physical features disappear for highly excited coherent states, where, however, the 
main quantum features are lost. 


2.6 Homodyne Detection 

This Subsection is devoted to the customary homodyne detector, which, despite it exits from the present 
phase estimation treatment, however it is the most relevant device in any interferometric setup. Actually, 
the homodyne detector belongs to the class of the zero-point measurement schemes, and thus is not a 
measurement of phase . The balanced homodyne scheme measures one quadrature of a field mode, which 
in turn is related to its phase difference with respect to the synchronous LO. Generally one is interested 
in the measure of the phase shift x of the signal state 

IVOx = exp(-txn)|^)o , (2.44) 

where, without loss of generality, the input state is assumed of the form 

oo 

IV-)o = £ c »l n ) , C„ > 0. (2.45) 

n=0 

The expectation value of the quadrature is given by 

oo 

(« <t>)x - X) V n + lCnC„+i COS (<f> - x) = (oo)o cos (<j> - x) . (2.46) 

n=0 

The quadrature a $ is proportional to the cosine of the phase with a proportionality “constant ” (do)o 
which can be evaluated from the knowledge of the fixed input state. Notice that, however, when the 
present scheme is regarded as a measure of the phase of the state |V») X itself the state-depending “constant 
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” is unknown, and it cannot be preventively measured without destroying the information on phase. 
In tlug case the only point which does not need any knowledge of the “constant” is the <f> — x = T /2 
point, namely just the maximum -derivative zero-current working point. Thus, it is not possible to detect 
a arbitrary phase difference, the sensitivity depending on the particular input signal, whereas a suitable 
feedback U needed to follow the working point. A convenient description of the homodyne 

detector in view of the above considerations is given in [17], where the zero-point (zero-field) probability 
distribution is revisited as a sort of a phase probability distribution. 

The phase sensitivity of the homodyne detector can be straightforwardly obtained from error propagation 
calculus using the relation 


6<!>=J(Kty 



(2.47) 


which is customary in the literature on interferometry. One can see that the $ - x = */ 2 working pomt 
minimizes sensitivity. 


3 Optimal states for phase measurements 


The design a phase measurement needs optimization of both the detection scheme and of the quantum 
state which carries the phase information. The former is the main task of quantum estimation theory, 
which lead s to an ideal scheme to be compared with the feasible ones. The latter, which is the main 
concern of this Section, depends on the detection scheme itself, and should account for the actual physical 
constraints, mainl y the total power impinged into the state. Therefore, the problem is that of optimizing 

the r.m.s sensitivity A<f> = f \J( A^ 2 ) for fixed average photon number, and depending on the particular 
detection scheme. 

In the following we consider, without loss of generality, a zero average phase state, with real coefficients on 
the number basis, namely 

i vo = 13 Cn i n ) Cn 6 R - ( 3>1 ) 

n—0 

The state optimization problem is to minimize a quantity of the form 

2 

<A^ 2 ) = + 2 13 A n<m c n c m , (3.2) 

6 n/m 


with the constraints 


13 C » = 1 » 13 « C n = <«) = » 


(3.3) 


n=0 


n=0 


The (real symmetric) matrix A = { A n m } depends on the detection scheme. In particular, for n ^ m one 
has (A ntn = 0) 


II 

;5 

(canonical) 

(3.4) 

r(“ + i) 

\/n!m! 

(DBH) . 

(3.5) 


A - v -J- 
n ' m (n - 

The method of Lagrange multipliers reduces the problem to that of minimizing the following expression 


^’({ c n}; A,/?|n) — -r- + 2 ^3 An ,m C n C m + 

n^m 


a(XX-i) +/3(f>;-n) (3.6) 
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with respect to {c n }, A and (3 being the Lagrange multipliers. The variational problem (3.6) is that of a 
quadratic symmetric form and is equivalent to the eigenvalue problem 

(M(/J) + AI) • c = 0 c •W (co, <*, ...) , (3.7) 

for symmetric matrix M = {Af nm } given by 

n,m — A nm -f 8 nm j3n . (3.8) 

Eq.(3.7) for the matrix (3.8) can be numerically solved upon suitable truncation of the Hilbert space 
H t . The absolute minimum corresponds to eigenvalue A = t 2 /3, and n 6 [0, dimW,/2] turns out to be 
a decreasing function of the running parameter /? € [0, 1]. Notice that one should consider only average 
values n <; dim7f,/2, such that the number distribution has vanishing tail at n = dim7f,, in order to avoid 
undesired numerical boundary effects. 



Figure 4: Phase and number probability distributions of an optimal states of n = 20 for both ideal and 
DBH detection 


3.1 Canonical Measurement 

For ideal measurement of the phase, the best phase states obtained through the above optimization 
procedure, lead to the simple power- law 


A 4> 




1.36 ±0.01 

fjl.00i0.01 ’ 


(3.9) 


in agreement with results of [18]. The proportionality constant actually increases very slowly as a function 
of n, and one has a variation of few percent for two decades of A <f>. Eq.(3.9) can be compared with the 
result of [18], and with the theoretical bound A <f> ~ l/(eh) [19] obtained by means of information-theory 
arguments. One should notice that essentially the same result can be obtained for large h (h > 10) using 
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squeezed states where the squeezing photon number is optimized as a function of the average total number. 
It turns out that the optimal states have only ~ 3.7% of squeezing photons (see Fig. 5). This result is 
quite different from the customary 50% optimal squeezing number (which also holds true for the homodyne 
sensitivity of the Mack-Zehnder interferometer [20]). 



<n> 


Figure 5: Opt imal squeezing photon number as a function of the average total number for both ideal and 
DBH detection. 


3.2 Double-homodyne 

As expected, an actual measurement of the phase does not achieve the ideal sensitivity (3.9). In the case 
of double-homodyne (or equivalently heterodyne) phase detection, the resulting power-law is 


AJ ( 1.00 ± 0 . 01 ) 

~ - 0 . 65 ± 0.01 ’ 


(3.10) 


which is obtained by numerically solving Eq.(3.7) for matrix M given in Eqs.(3.8) and (3.5). In Fig. 4 the 
optimized states for both canonical and DBH detection are compared for an equal fixed average photon 
number n = 20. One ran see that the number and phase probability distributions are qualitatively similar, 
however the DBH optimum states are slightly sharper in the number distribution and larger in the phase 
one. Also the best DBH states are essentially indistinguishable from squeezed states which are optimized 
in the squeezing photon number as a function of the average total number (see Fig. 5). In this case, only 
less than ~ 2% of squeezing photons turns out to be optimal. 
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<n> 


Figure 6: Uncertainty product versus the average photon number for optimal states. 

3.3 Heisenberg Uncertainty Product for phase quadrature 
The customary Heisenberg uncertainty relation 

AAAB>±\Tt{p[A,B]}\ (3.11) 

refers to the situation in which the quantum system is prepared in a state with fixed uncertainty say A<4 
and the other observable B is measured. For the case of a joint A-B measurement, however, a generalized 
uncertainty relation hold, where the 1/2 factor on the right side of Eq.(3.11) is dropped, corresponding to 
an added noise of 3 dB [21]. For the phase-quadratures one has the commutation relation 


M = ~l°><o|, 

(3.12) 

corresponding to the joint-measurement uncertainty relation 


AcA a >i|(V>|0)| 2 . 

(3.13) 

In Eqs.(3.11) and (3.13) the uncertainties are defined in the usual way, namely AO 2 
other hand, in the POM approach the correct uncertainty (namely the measured 

AO 2 = {O 2 ) — (0) 2 , where 0 2 ^ 0 2 is defined as in Eq.(2.4). In general, by means 
one obtains 

= ( 6 2 ) - (O) 2 . On the 
quantity) is defined as 
of Schwartz inequality, 

(o 2 ) - (6) 2 > (O 2 ) - <0) 2 . 

(3.14) 
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Ideal detection 

Double-homodyne 

Homodyne 

Ideal Sine detection 


BBfflSiffli 


A* (Eq.(2.47)) 

■BBBBi 

Coherent 

SUtco 



ift-i 

ifi-i 

Phan* Coherent 
States 
(see [5]) 

(ln4)iir* 



ifi-i 

Optimised 

Squeezed 

States 

(i.je±o.oi)ft -,00:t ° l 

»,,/* ~ 4% 

(1.00 ±0.01)ft~° 851: 01 
~ 2% 

ift - 1 

n J(? /n = 50% 

(1.36 ±0.01)ft _1 001: 01 
litq/li ~ 5% 


Table 1: Asymptotic sensitivities versus the average photon number (results from numerical calculations 
are given with error estimation) 

Therefore, one should not expect that the optimum states for phase detection achieve the minimum 
uncertainty product (3.13), even though phase detection corresponds to a joint sine-cosine measurement. 
In Fig. 6, the minimum uncertainty product in Eq.(3.13) is compared with the actual uncertainty product 
of the optimal states for canonical measurement. One can see that the minimum uncertainties are 
never achieved, despite canonical detection is ideal. In the same figure the uncertainty product of the 
optimal states for DBH detection is also reported for comparison, showing additional noise due to nonideal 
measurement. 

3.4 Homodyne detection 

For homodyne detection the well known sensitivity A <f> = ^n” 1 is achieved only near the zero-current 
working point. As a consequence, it happens that it is better than the sensitivity achieved by the true 
phase measurement, either in the double-homodyne or the ideal case itself. This is due to the fact that the 
measurement of the field- quadrature near the zero-current working point partially underestimates the tails 
of the phase distribution at <f> = ±tt. The latter are enhanced by large squeezing, and thus one also finds 
that the optimal number of squeezing photons is only a few percent of n for the true phase measurements, 
whereas it is 50% for single homodyning. It is interesting to notice that the measurement of a single 
phase-quadrature, say sin^, also exhibits a small optimal squeezing fraction (~ 5%), as, in some sense, it 
is more faithful observable than the field quadrature. In Tablel the above results are reported along with 
the phase sensitivity for other quantum state and different detection schemes in the limit of large average 
photon numbers n. 

4 Conclusion 

We have analyzed ideal and actual detection schemes for the quantum phase. We have also considered the 
homodyne detectors as it can be used as zero-phase measurement apparatus. The working conditions for a 
phase detectors have been discussed showing that in an actual measurement the phase shift corresponds to 
the polar angle between two real output photocurrents. We have analyzed in detail the double-homodyne 
scheme of [11], giving the POM of the apparatus. The equivalence of this schemes with heterodyne one in 
detecting the phase it has also been shown. 

A critical revision of various adopted definitions of sensitivity has been reported, we have concluded that 
the usual r.m.s. noise is the right quantity to be considered. 
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The sensitivity of all the detection schemes has been optimized at fixed energy with respect to input 
quantum states of radiation. We have shown that the sensitivity versus the average photon number n is 
bounded by the ideal limit A^ ~ n -1 , whereas for double-homodyne detection the bound is A <f>o ~ n“ 2 / 3 , 
in between the shot noise level A <f> ~ fi *^ 2 and the ideal bound. The optimal states achieving the best 
sensitivity for fixed energy have been numerically obtained, and we have shown that they are very close 
to coherent states wealdy squeezed either for ideal detection or for double-homodyne. This result, which 
is in contrast with the 50% of optimal squeezing photons for single- homodyning, is due to the sensitivity 
of the double-homodyne detection to the whole phase distribution including tails, which are enhanced by 
increasing squeezing. 

Appendix: Evaluation of Double-Homodyne POM 
The phase distribution is the marginal probability integrated over the modulus p , namely 


p {4) = J q pdp P\ (p cos <f>)P 2 (p sin $) . (A.l) 

Using Eq. (2.21) one has 

poo dti dv 

= / pdp / — tl{pp ® ps e'* i ( ri_pcos ^) +, * , ( J 2 -p 8 in^)} (A.2) 

Jo J-r\x\ 2t |,| 2x J K ' 

ps being the density matrix of the mode a (the system) and 

PP = |0)(0| ® |0)(0| ® \z){z\ (A.3) 

the density matrix of the probe. From Eqs. (2.6) and (2.7) one can see that the POM is obtained upon 
tracing over the probe Hilbert space “Hp, thus obtaining the operator which acts on the system space 7is 
only 

POO pk\z\ ptt\z\ 

dpo{<t>) = d<j> pdp ^ — tT P {p P ® i s e'^-pcos^+i^ij-psin^, (A.4) 

JO J — ir|z| J — i r|z| 

Using the coherent state resolution of identity one has 

dpD{<f>) = d<t> f pdp f 11 ^ ^ / — f ~e~ ip ^ cos,t>+, ' sia ^\u)R{w \ , (A.5) 

Jo J-k\z\ 2t J- t |,| 2ir Jc ir Jc * v ' 

where R is the matrix element 

R= (u,z t 0 t 0\UV^ + ^U\0,0 t z,w). (A.6) 

In Eq. (A.6) U denotes the unitary evolution operator of the detector, which acts on the state (A.6) as 
follows 

U\0,0,z,w) = \^(z + w),^{z - w),^(iz + w),^(iz - w )) . (A.7) 

The explicit expression of the matrix element R is given by 

li* = H^-^H ! -iwHiw 2 (co S i+cos^) + izfi(ism^+i,m^) 

+ ^(jsm^-i.ini)+l»fi(co.^ + c<> s ^) . (A.8) 
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Taking the «trong LO limit \z\ -* oo and introducing the complex variable a = |(v + s#*)e ,ar *W one get» 

OO i 

R = exp f--|«| 2 - ;;M 2 - l«| 2 + au- tua + tim = — {u|a t |a)(“°l®l u ') • ( A - 9 ) 

12 2 J p=0 P' 

Substituting Eq.(A.9) into Eq.(A.5) leads to 


r f" -e~^ 


costfr+v sin 


2x 


*>V± [ — |«)(«|d-*|c){-a|a» / ^|u,){a.| 

^ 0 pl Jc * JC * 


«p 


dM*) = 

. d*f'~i* [2 r p dpf — « p 0 ^»)| o >( o |^** «*(/***-») 

"p!2x [ Vo /C » J 

Using the coherent resolution of the identity and integrating over p one obtains 

dun(6) JV- f , 

2x “.?! ^“V n!m! 

p— 0 ^ m,n=0 

where T(*) is Eulero’s Gamma function. The normal ordered representation of the vacuum state 

lim £ ^T^a pt “ P = l°>{°l. 

e ^ P =o P ! 


(A.10) 


(A.11) 


leads to 


dp D {<f,) = ^ f; (_.)n-m e ,(n-mM r ( Ii f^+ ]} e *f" a"* t |0>(0|« n e'f" . 

m,n=0 n.m. 

Prom Eq.(A.13) one obtains the POM of the detector in form of a double series 

4 M « = g E '"><”*1 • 

n,m v 

Alternatively, using the T-function integral representation one can write 

dp D {4>) = ^ [°°pdp e“'^e pe * at |0){0|e pe ^ [°°pdp \pe»)(pe *\ . 

T ;o » 0 


(A.12) 


(A.13) 


(A.14) 


(A.15) 


146 



References 

[1] G. M. D’Ariano, and M. G. A. Paris, submitted to Phys. Rev. A 

[2] P. A. M. Dirac. Proc. R. Soc. London A 114 243 (1927) 

[3] D. T. Pegg, S. M. Barnett Phys. Rev. A39, 1665 (1989) 

[4] C. W. Helstrom Quantum Detection and Estimation Theory (Academic Press, New York, 1976) 

[5] J. H. Shapiro, S. R. Shepard Phys. Rev. A43, 3795 (1991), J. H. Shapiro, S. R. Shepard, N. C. Wong 

Phys. Rev. Lett. 62, 2377 (1991), J. H. Shapiro in The Workshop on Squeezed States and Uncertainty 
Relations ed. by D.Han et al. (NASA Conf. Public. 3135, Washington DC, 1992), p.107 

[6] A. S. Lane, S. Braunstein, C. M. Caves Phys. Rev. A47, 1667 (1993) 

[7] L. Susskind and J. Glogower, Physics 1 (1964) 49. 

[8] M. Ozawa, ‘Operator algebras and nonstandard analysis’ in Current Topics in Operator Algebras ed. 

by H. Araki et. al., World Scientific, (Singapore 1991) p.52 

[9] M. Ban J. Opt. Soc. Am. B9, 1189 (1992) 

[10] H. P. Yuen Phys. Lett. 91A, 101 (1982) 

[11] J. W. Noh, A. Fougeres, L. Mandel Phys. Rev. Lett. 67, 1426 (1991); Phys. Rev. A45, 424 (1992); 
Phys. Rev. A46, 2840 (1992) 

[12] Z. Hradil Quantum Opt. 4, 93 (1992) 

[13] A. Bandilla, H. Paul, H. H. Ritze Quantum Opt. 3, 267 (1991) 

[14] Z. Hradil, J. H. Shapiro Quantum Opt. 4, 31 (1992) 

[15] P. Carruthers, M. M. Nieto Rev. Mod. Phys. 40, 411 (1968) 

[16] M. J. W Hall, I. G. Fuss Quantum Opt. 3, 147 (1991) 

[17] W. Vogel, W. Schleich Phys. Rev. A44, 7642 (1991) 

[18] G. S. Summy, D. T. Pegg, Opt. Commun.77,75 (1990) 

[19] H. P. Yuen, in The Workshop on Squeezed States and Uncertainty Relations ed. by D.Han et al. 
(NASA Conf. Public. 3135, Washington DC, 1992), p.13 

[20] R. S. Bondurant and J. H. Shapiro, Phys. Rev. A30 2548 (1984) 

[21] E. Arthurs, M. S. Goodman Phys. Rev. Lett. 60, 2447 (1988) 


147 




